Universal impurity-induced bound state in topological superfluids 
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We predict a universal mid-gap bound state in topological superfluids, induced by either non- 
magnetic or magnetic impurities in the strong scattering limit. This universal state is similar to the 
lowest-energy Caroli-de Gennes-Martricon bound state in a vortex core, but is bound to localized 
impurities. We argue that the observation of such a universal bound state can be a clear signature 
for identifying topological superfluids. We theoretically examine our argument for a spin-orbit 
coupled ultracold atomic Fermi gas trapped in a two-dimensional harmonic potential, by performing 
extensive self-consistent calculations within the mean-fleld Bogoliubov-de Gennes theory. A realistic 
scenario for observing universal bound state in ultracold *"K atoms is proposed. 
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Topological superfluids are of great interest [1]. They 
are promising candidates that host Majorana fermions 
[2II, which lie at the heart of topological quantum infor- 
mation and computation, due to their exotic non-Abelian 
exchange statistics [^-Q. To date, there are a number 
of proposals for practical realizations of topological su- 
perfluids, including p + ip superconductors [6, 7], sur- 
faces of three-dimensional (3D) topological insulators [8- 
llOl or one-dimensional (ID) spin-orbit coupled nanowires 
[ill . [1^ in proximity to an s-wave superconductor, and 
two-dimensional (2D) (isHiit or ID [IMl spin-orbit 
coupled atomic Fermi gases near Feshbach resonances. 
All these proposals are appealing and are to be examined 
experimentally. In fact, recent experimental results on 
tunneling spectroscopy of semiconductor InSb nanowires 
in a magnetic field placed in contact with a supercon- 
ducting electrode [20| may already suggest the existence 
of topological superfluids and Majorana fermions. How- 
ever, unambiguous characterizations of topological prop- 
erties of the nanowires are still missing. 

In this Letter, we propose that a universal mid-gap 
bound state, induced by strong non-magnetic or mag- 
netic impurity scattering, could provide a clear signa- 
ture for the existence of topological superfluids. In solid 
state, impurities are widely known to serve as an im- 
portant local probe that characterizes the quantum state 
of hosting systems (2l| . Individual impurities have used 
to determine the superconducting pairing symmetry of 
unconventional non- s-wave superconductors [i^l and to 
demonstrate Friedel oscillations on Be(OOOl) surface [l^. 
In strongly-correlated many-body systems, they may be 
employed to pin one of the competing orders j24|. Here, 
unique to topological superfluids, we predict that a sin- 
gle impurity with sufficiently strong scattering strength 
can create a universal mid-gap state bound to the im- 



purity. It resembles the lowest-energy Caroli-de Gennes- 
Martricon (CdGM) bound state inside a vortex core [2^. 
For small order parameters, where the bound state en- 
ergy E is nearly zero, the wave- function of the universal 
bound state is found to closely follow the symmetry of 
that of Majorana fermions fl^. 

In our work, the emergence of universal impurity- 
induced bound state is examined theoretically in an in- 
teracting spin-orbit coupled ultracold atomic Fermi gas 
in 2D harmonic traps [16]. We perform numerically ex- 
tensive self-consistent calculations by using fully micro- 
scopic Bogoliubov-de Gennes (BdG) theory, to explore 
the details of the universal bound state. This specific 
choice of topological superfluids is motivated by the re- 
cent realization of spin-orbit coupling in atomic Fermi 
gases of [2^ and ®Li atoms Benefited from the 
high controllability in interaction, geometry and purity 
in cold-atom experiments, 2D spin-orbit coupled atomic 
Fermi gases are arguably the best candidate for observing 
the predicted universal bound state. Our results, how- 
ever, should be applicable as well to various topological 
superfluids that are believed to exist in solid state. We 
propose a realistic scenario of creating universal bound 
state in ^'^K atoms and discuss briefly the relevance of 
our results to other solid state systems. 

Mean-field BdG equation. — To start, we consider a 
trapped 2D atomic Fermi gas with a Rashba-type spin- 
orbit coupling and a Zeeman field h, which is believed to 
be a topological superfluid when the Zeeman fleld exceeds 
a threshold he [3| ■ The model Hamiltonian of the system 
is given by, 'H = / dr['Ho(r) + 'H/(r) + ■Himp(r)], where 

Ho(r) = ^i^5(r)V'a + [i^lVso{r)^i + H.c] (1) 

is the single-particle Hamiltonian density in the presence 
of Rashba spin-orbit coupling Vso(r) = —i\{dy -\- idx), 
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■H/(r) = Uoip^{r)^jj]^{r)ipi{r)ipf{r) represents the inter- 
action, and ?^imp(r) = E^=t,i V'i^^imp('')V'<T describes 
the potential scattering due to the impurity. Here, ip'!^ ^ 
are respectively the creation field operators for the spin- 
up and spin-down atoms and, T-L^{r) = — ?i^V^/(2M) + 
MuPr"^ /2 — /i — haz is the single-particle Hamiltonian 
in a 2D harmonic trapping potential MitP'r^ jl. in ref- 
erence to the chemical potential /i. We have used the 
standard s-wave contact interaction between atoms with 
opposite spins, whose strength [/q is to be regularized 
by the binding energy of the two-body bound state Ea 
liL 28| . For computational simplicity, we place an impu- 



rity at origin and consider either a delta-like scattering 
potential, V5J^p(r) = l^^pj(r), or a gaussian-shape poten- 
tial with width d, Y-^^Jyv) = [\^^p/(7rd2)]exp[-rVd2]. 
In the case of magnetic impurity, we take the potential 
strength T^^p = — V^^^ = — Mmp; while for non-magnetic 

impurity, V^^^ = V^j^^p = — Mmp- We have checked both 
positive and negative values of Vlmp and have observed 
very similar results at large |Viinp|. Hereafter, we focus 
on the case with Vlmp > 0. 

We solve the low-energy fermionic quasiparticles of 
the model Hamiltonian by using the standard mean-field 
BdG approach, 'H.^ag^^] (r) = -E',,^'^ (r), where 
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is the BdG Hamiltonian, ^'r; (r) = [u-f-^, it^,,, w^,,]"^ 
and Erf are the Nambu spinor wave-functions and en- 
ergies for quasiparticles, respectively. Within mean- 
field, the order parameter takes the form A(r) = 
-(C^o/2)E^[wt'7«|^/(-E^r,) and, is to be 

solved self-consistently together with the atomic densi- 
ties, n,(r) = (1/2) + 1^-^,1 

Here / (x) = 1/ (e^/''^^ -t- l) is the Fermi distribu- 
tion function at temperature T. The chemical po- 
tential /X, implicit in 'H^(r), can be determined by 
the total number of atoms iV using the number equa- 
tion / (ir[n-t- (r) + (r)] = N . As the impurity is 
placed at origin r = 0, the BdG Hamiltonian pre- 
serves rotational symmetry. Therefore, we take A(r) — 
A(r) and decouple the BdG equation into different an- 
gular momentum channels indexed by an integer m, 
with which the quasiparticle wave functions become, 
[ut„(r),u^^(r)e*'^,^;t^(r)e*^,^;4^(r)]e™'^/v^. By ex- 
panding Uari (f) and v^r; (f) in the basis of 2D harmonic 
oscillators, the solution of BdG equation converts to a 
matrix diagonalization problem. Numerically we have 
to truncate the summation over energy levels 77. This 
is done by introducing a high energy cut-off Ec, above 
which a local density approximation is used for high-lying 
wave- functions [29] . We have checked that such a hybrid 
procedure is numerically very efficient. 

For the results presented here, we have solved 
self-consistently the BdG equation for a cloud with 
N = 400 atoms at zero temperature. In 2D har- 
monic traps, it is convenient to use the Fermi ra- 
dius rp = {4:Ny/^y/h/{Muj) and Fermi energy Ep — 
fi^kp/ {2M) — \/Nhjjj as the units for length and en- 



ergy, respectively. The strength of impurity scattering 
potential V^^-^^ will be measured in units of rpEp. We 
have taken an interaction parameter Ea = Q.2Ep and 
a spin-orbit coupling strength Xkp/Ep — 1. With these 
parameters, the whole Fermi cloud becomes a topological 
superfiuid when the Zeeman field is larger than a thresh- 
old he ~ O.blEp [1^. Let us first consider the localized 
impurities with a delta- like scattering potential Vj^pS{r). 

Emergence of universal impurity bound state. — Ac- 
cording to Anderson's theorem [30], a conventional s- 
wave superfiuid can barely be affected by non-magnetic 
impurities. In contrast, magnetic impurities can break 
time-reversal symmetry of superfiuid and act as pair 
breakers, leading to the appearance of mid-gap state - 
the so-called Yu-Shiba state - which is bound to local- 



ized impurities inside the pairing gap 3l|, |32| . The en- 
ergy of such a mid-gap bound state is determined by 
the strength of the impurity scattering potential VSmp. 
As Vimp increases, the Yu-Shiba state moves from the 
upper gap edge to the lower gap edge for the spin-up 
atoms and moves oppositely for the spin-down atoms. In 
the presence of Rashba spin-orbit coupling, we have con- 
firmed numerically that the above statements continue 
to hold, even under a Zeeman field, if the Fermi cloud 
is not a topological superfiuid. For a typical parameter 
h — 0.2Ep, with increasing the strength of magnetic im- 
purity, we find that the position of the Yu-Shiba state 
moves very quickly from one gap edge to the other. 

In contrast, once the Zeeman field is beyond the thresh- 
old he so that the whole Fermi cloud becomes a topo- 
logical superfiuid, we observe entirely different behavior, 
as revealed in Fig. 1. For non-magnetic impurities, an 
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FIG. 2: (color online) Gap parameter as a function of im- 
purity strength Vimp (in units of r^Ep), for a non-magnetic 
impurity (a) and for a magnetic impurity. In the limit of 
strong impurity scattering, the gap parameter has the same 
spatial distribution, no matter the impurity is non-magnetic 
or magnetic. 



FIG. 1: (color online) Bound states induced by a non- 
magnetic delta-like impurity (a) and by a magnetic delta- 
like impurity (b), Vi"-^p{r) — Vi"-^-^pS{r), in a topological su- 
perfluid at /i = Q.lEp, as shown by the peaks in the to- 
tal local density of state (LDOS) p{r,E) at kpr — 2. Here, 
P{r, E) = P<.(r, E) and p^r, E) = (1/2) EJh<.^l' 5{E - 
Er,) 4- \var]\^ 5{E + Er,)]. The dashed and dash-dotted lines 
highlight the resonance peak position or the energy of bound 
states. From bottom to top, the impurity strength increases 
from Vimp = to Viu,p = 0.011r|-£;F, in steps of O.OOlrf-^F. 
The curves are offset for clarity, except for the lowest curve 
at Vimp = 0. (c) The energy of bound states as a function of 
the impurity strength, in units of the gap parameter at the 
trap center in the absence of impurity, Ao — G.2iQ1Ef- 

unexpected bound state appears from one gap edge as 
the impurity strength is larger than a critical strength 
Mmp ^ Q.QQAr^Ep- As l^mp increases, the bound state 
moves towards, but never reaches zero energy. In fact, 
its energy saturates quickly to i? ~ Q.llEp — Ag/£'i?, 
where Aq ~ 0.307i?F is the gap parameter at the trap 
center in the absence of impurity. For magnetic impuri- 
ties, the dependence of the position of the Yu-Shiba state 
on the impurity strength is also strongly modified: at 
large impurity scattering, the Yu-Shiba state now moves 
to ~ Aq/E^f, nearly at the same energy as the new 
bound state induced by strong non-magnetic impurities. 
This coincidence in the energy of bound states clearly in- 
dicates that in topological superfluids a universal bound 
state emerges in the limit of strong impurity scattering. 

Origin of the universal state. — The appearance of 
bound states implies that the gap parameter would be 
strongly depleted close to the impurity. In Fig. 2, we 
examine the spatial profile of order parameter near the 
impurity. For a weak non-magnetic impurity, as shown 
in Fig. 2(a), the gap parameter is already strongly modi- 



fied at T^mp ^ 0.004r|,i?i?. Being regarded as a scattering 
potential for Bogoliubov quasiparticles [2^, the gap pa- 
rameter hence starts to accommodate a bound state. For 
a weak magnetic impurity (Fig. 2(b)), the pair-breaking 
effect is always significant enough to induce a Yu-Shiba 
bound state, as anticipated. In the strong scattering 
limit, it is remarkable that the gap parameter acquires 
a universal spatial profile, despite the type and strength 
of impurities. It is fully depleted at the impurity site 
and has a very similar distribution as the gap parameter 
inside a vortex core. Therefore, we anticipate that the 
observed universal bound state would resemble the well- 
known CdGM vortex-core bound states [2^ . Indeed, the 
energy of the universal impurity state, E ~ Ag/Ep, is at 
the same order as that of CdGM bound states. 

Now, the formation of the universal bound state can 
be easily understood from its analogy with the CdGM 
vortex-core state. As the gap parameter is fully sup- 
pressed at the impurity site, we have a local point de- 
fect (i.e., vacuum) that is topologically trivial. Due to 
the topological nature of the Fermi cloud away from 
the impurity, there would be an interface between the 
non-topological and topological com pon ents, which can 
host a gapless Majorana edge state [33j. The observed 
universal impurity state is precisely such a Majorana 
edge mode. However, its energy is not ex actly zero due 
to the finite confinement of the system [3J|. As de- 
rived analytically by Stone and Roy [s^ (see also Ref. 
[3^). the dispersion relation of edge states in topologi- 
cal superfluids with a confinement length ^ is given by 
E{m) = —(to + 1/2)Ao/(A:fC)- By assuming a charac- 
teristic length ^ ^ hvp/Ao for the gap parameter distri- 
bution (25| . where vp the Fermi velocity, we estimate 
that E ~ Aq/Ep, in good agreement with the observed 
energy of the universal bound state. 
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FIG. 3: (color online) The wave-function of the universal 
bound state, induced by either non-magnetic or magnetic im- 
purity in the strong scattering limit. The inset shows the lin- 
ear contour of LDOS for spin-up and spin-down atoms near 
the impurity site. 
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FIG. 4: Loss of the universal bound state for an extended 
impurity. Here, we take a gaussian-shape scattering potential, 
^4mp(r) = [\/iSip/(7rd^)lexp[-rVd^], with width kpd = 0.5. 
From bottom to top, the impurity strength increases from 
Vimp = to Viu,p = QmrlEp, in steps of 0.002r|,£;F. Other 
parameters are the same as in Fig. 1. 



In Fig. 3, we examine the wave-function of the univer- 
sal bound state. Indeed, it satisfies approximately the 
symmetry Mg. (r) = u* (r) , which should be obeyed by 
zero-energy Majorana fermions. In the inset, we present 
the LDOS close to the impurity site. The universal bound 
state is clearly visible within the gap. Experimentally, 
the LDOS may be measured through spatially resolved 
radio-frequency (rf) spectroscopy [36], which provides a 
cold-atom analog of the widely used scanning tunneling 
microscope in solid state (s^l- The wave- function of the 
universal bound state can therefore be determined from 
the real-space structure of LDOS within the gap. 

Loss of university. — The universality of the impurity- 
induced bound state can be lost if the impurity scattering 
has a finite width. In this case, a hole will be created in 



the strong impurity scattering limit, instead of a point 
defect. Therefore, there are a series of edge states. The 
wave-function and energy of these edge states would de- 
pend critically on the shape and strength of the impurity 
potential. In Fig. 4, we show the bound states induced by 
a non-magnetic (a) and a magnetic (b) gaussian impurity, 
with a finite width kpd = 0.5. It is readily seen that with 
increasing the impurity strength the bound state never 
approaches to a universal limit. We have checked that 
for larger widths, the LDOS becomes very complicated, 
as more and more bound states appear. 

Experimental proposal. — We now show that ultra- 
cold Fermi gases of atoms is a potential candidate 
for observing the predicted universal impurity-induced 
bound state. A 3D spin-orbit coupled Fermi gas was 
recently realized at Shanxi university (26| . By loading 
a pancake-like optical trap V{v,z) = M[a;^r^ -I- a;^z^]/2 
with trapping frequencies a;^ ^ a; (38| or using a deep 
ID optical lattice [HI, a 2D topological superfluid with 
number of atoms N ~ 1000 and size rp ^ lOO/im may 
be prepared at the temperature about lOnK . It is conve- 
nient to create the delta-like impurity potential by using 
a dimple laser beam that has a sufficiently narrow beam 
width d < 1/iTO [4^, so that kpd <^ 1. By suitably 
tuning its frequency, the scattering potential caused by 
the laser beam can be attractive or repulsive for different 
spins. Thus, both non-magnetic and magnetic impurities 
can be simulated. The resulting universal bound state 
may be visualized by using the standard tool of spatially 
resolved rf-spectroscopy. All the techniques required to 
observe the predicted universal state are therefore within 
the reach of current experiments. 

Application to other solid state systems. — Our results 
are apparently applicable to the triplet superconductor 
Sr2Ru04. For the possible ID topological superconduc- 
tor reported recently in InSb nanowires [20|, a strong 
impurity potential would split the ID topological super- 
conductor into two. Therefore, at the impurity site we 
anticipate two universal bound states, with precise zero- 
energy. The observation of such a pair of zero-energy 
Majorana fermions is an unambiguous identification of 
the topological nature of InSb nanowires. 

Conclusion. — We have investigated the non-magnetic 
and magnetic impurity scattering in an atomic topolog- 
ical superfiuid and have predicted the existence of uni- 
versal bound state for strong impurity scatterings. The 
observation of such a universal bound state - via spatially 
resolved radio- frequency spectroscopy - is a smoking-gun 
proof of atomic topological superfluidity. Our prediction 
seems within experimental reach and opens the way to 
unambiguously characterizing the topological properties 
of other solid-state systems, such as the unconventional 
superconductor Sr2Ru04 and ID topological supercon- 
ductor of InSb nanowires. 

Note added. — After completing this work, we were 
aware a related non-self-consistent T-matrix calculation 
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in ID topological superconductors, which predicted a 
bound state induced by non- magnetic impurities [41]. 
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